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Abstract 

In this article we describe a model of the universe consisting of a mixture of the or- 
dinary matter and a so-called cosmic quaternionic field. The basic idea here consists 
in an attempt to interpret A as the energy density of the quaternionic field whose 
source is any form of energy including the proper energy density of this field. We set 
the energy density of this field to A and show that the ratio of ordinary dark matter 
energy density assigned to A is constant during the cosmic evolution. We investigate 
the interaction of the quaternionic field with the ordinary dark matter and show 
that this field exerts a force on the moving dark matter which might possible create 
the dark matter in the early universe. Such determined A fulfils the requirements 
asked from the dark energy. In this model of the universe, the cosmical constant, the 
fine-tuning and the age problems might be solved. Finally, we sketch the evolution 
of the universe with the cosmic quaternionic field and show that the energy density 
of the cosmic quaternionic field might be a possible candidate for the dark energy. 



1 Introduction 

According to the astronomical observation, the evidence continues to mount that the 
expansion of the universe is accelerating rather than slowing down. New observation 
suggests a universe that is leigh- weight, is accelerating, and is flat [14] [12] [13]. To 
induce cosmic acceleration it is necessary to consider some components, whose equations 
of state are different from baryons, neutrinos, dark matter, or radiation considered in 
the standard cosmology. 
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As is well-known, in cosmology a new kind of energy is considered called quintessence 
("dark energy"). Quintessence represents a dynamical form of energy with negative 
pressure [6]. The quintessence is supposed to obey an equation of state of the form 



p Q c 2 = WQQQ, 



-Kwq < 0. 



(1) 



For the vacuum energy (static cosmological constant), it holds wq = — 1 and wq = 0. 

An adequate cosmological theory, conform with recent observation, should give an- 
swers to the following problems [39] [40]: 

(i) The cosmological constant problem. The 'A-problem' can be expressed as discrepancies 
between the negligible value of A for the present universe and the value 10 50 times larger 
expected by Glashow-Salam- Weinberg model [1] or by GUT [2] where it should be 10 107 
times larger. 

The fine-tuning problem. It is a puzzle why the densities of dark matter and dark en- 
ergy are nearly equal today when they scale so differently during the expansion of the 
universe. Assuming that the vacuum energy density is constant over time and the mat- 
ter density decreases as the universe expands it appears that their ratio must be set to 
immense small value (« 10~ 120 ) in the early universe in order for the two densities to 
nearly coincide today, some billions years later. 

(hi) The age problem. This problem expresses the discrepancy connected, on the one 
side, with the hight estimates of the Hubble parameter and with the age of globular 
clusters on the other side. The fact that the age of the universe is smaller than the age 
of globular clusters is unacceptable. 

(iv) The flatness problem. Inflation predicts a spatially flat universe. According to Ein- 
stein's theory, the mean energy density determines the spatial curvature of the universe. 
For a flat universe, it must be equal to the critical energy. The observed energy density 
is about one-third of critical density. The discrepancy between the value of the observed 
energy density and the critical energy is called the flatness problem. 

(v) Problem of the particle creation. In variable A models the creation of particles gen- 
erally takes place. The question what is the mechanism for this process represents the 
problem of the particle creation. 

It is well-known that the Einstein field equations with a non-zero A can be rearranged 
so that their right-hand sides consist of two terms: the stress-energy tensor of the ordinary 
matter and an additional tensor 



A is identified with vacuum energy because this quantity satisfies the requirements asked 
from A, i.e. (i) it should have the dimension of energy density, and (ii) it should be 




(2) 
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invariant under Lorentz transformation. The second property is not satisfied for arbi- 
trary systems, e.g. material systems and radiation. Gliner [4] has shown that the energy 
density of vacuum represents a scalar function of the four-dimensional space-time coor- 
dinates so that it satisfies both above requirements. This is why A is identified with the 
vacuum energy. 

From what has been said so far it follows that the following properties are required from 
the vacuum energy density: (i) It should be intrinsically relativistic quantity having the 
dimension of the energy density, (ii) It should be smoothly distributed throughout the 
universe, (iii) It should cause the speedup of the universe, (iv) It should balances the 
total mean energy density to £1 = 1. 

In the next Sections we will describe a model of the universe consisting of a mixture 
of the ordinary matter and a so-called cosmic quaternionic field. The basic idea in this 
article consists in an attempt to interpret A as the energy density of the quaternionic field 
whose source is any form of energy including the proper energy density of this field The 
article is organized as follows. In Section 2 we describe the proposed quaternionic field. 
In Section 3 we set the energy density of this field to A and show that the ratio of ordinary 
dark matter energy density assigned to A is constant during the cosmic evolution. In 
Section 4 we show that this field exerts force on moving dark matter. In Sections 5 we 
describe the possible mechanism of the particle creation. We sketch the evolution of the 
universe with the cosmic quaternionic field. Finally, we show that such defined A fulfils 
the requirements asked from the dark energy and that in this model of the universe the 
above problems might be solved. 

2 The cosmic quaternionic field 

In a very recent article [11], A has been interpreted as the field energy of a quaternionic 
field (called <I>-field, for short) [5] [10] [9] (see also the Appendix). The field equation of 
the $-field can be written in the following form 

diFij = Jj, (3) 

where Jj = kpvi ( Vi being the 4- velocity, p matter density) is the current of the ordinary 
matter and Jo = k(e se if + p)v$. e se if is the energy density of the <£-field and Fy is the 
field tensor defined as 
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Fij has only diagonal components Fa = <I> i = 1,2,3, Fa = — $ i = 0, and Fij = 
i 7^ j. These components are transformed under Lorentz transformation as follows 
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[3] 

F33 = F33 = $ F M = ^2 (^44 + /^ll) = = F44. 

The field variable <D has the dimension of field strength and its square has the dimension 
of energy density. Two scalars can be formed from i 7 ^-: its trace F\ and F^F 13 the later 
has the dimension of the energy density. In the differential form the field equations (3) 
are 

V<Z> = kJ = kpVi, i = l,2,3 (4) 

and 

1 <9$ 

~--<ft = k o{e S elf + P)vi, i = (5) 
where e se i / is the energy density of the $-field given as [5] 

tself = $ 2 - 

From the Newton gravitation law it follows that the gravitational "charge" of a point 
mass rriQ is \/~GmQ [37] [23]. Accordingly, we set for the coupling constants k and ko 

k = and k$ = VG. 

c 

These equations are first-order differential equations whose solution can be found 
given the source terms. Assuming the spacial homogeneity of the $-field and the absence 
of any ordinary matter, i.e. J\ = J2 = J3 = Jo = 0. , the field variable <D becomes 
independent of spatial coordinates. Here, the only source of the $-field is its own energy 
density, i.e. e se if = <J> 2 . Therefore, Eqs.(4) and (5) become 

V$ = (6) 

1 d$ _ Vg$ 2 _ \/G$ 2 

c dt c 2 c 2 
The solution differential equation (7) has a simple form 



(7) 



where to is the integration constant given by the boundary condition. The energy density 
e se if is approximately equal to the observed value of the contemporary cosmological 
constant. 
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3 The dark energy modeled by a time-dependent cosmo- 
logical constant 



The theory of the time-dependent cosmological constant in the Friedmann model is well 
established ( see, e.g., [16]). The time-dependent cosmological models in the framework 
of scalar field theory were first discussed by P. J. E. Peeble and B.Ratra [20], B.Ratra 
and P. J. E. Peebles ([17]) and M. Ozer and M. O. Taha [18]. A number of authors 
set phenomenologically A oc 1/t 2 [22-29] (for a review see [21]). Generally, A contains 
in its definition the gravitation constant G and velocity of light c. From the purely 
phenomenological point of view the simplest expression for A oc 1/t 2 , having the right 
dimension, and containing G and c is 



A 



K 2 C 2 



8irGt 2 ' 

where ft is a dimensionless constant. Because A is set equal to the field energy density 
of the cosmic quaternionic field we have (to = c = 1) 

1 1 



A -5 



VGtVGt 



$2 1 



Accordingly, we have 



<J>2 i i 
A = — = — -—5- and A = (10) 
8tt 8vrGt 2 t 2 V ' 



The gravitational field equations with a cosmological constant A and the energy conser- 
vation law are (k=0) 

* = + A, H=* A=^ (11, 

R 4irG , , . 

^ = -^-(p + 3p + 2A) (12) 

and 

p + 3^(p + p) = -L (13) 
Suppose we have a perfect-gas equation of state 

p = ap (14) 

and suppose that the deceleration parameter is constant. If the evolution of the scale 
factor is given in form R oc t n then q = — (n — l)/n, therefore, we set 

1 - r*~\ 

q = =— = 1. (15) 

R 2 n y ' 
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If we suppose the time dependence p and A in form 

A B 
p = and A = B = const. and A = const., (16) 

respectively, then, inserting (16) into (12), (13) and (14), gives the following relation 
between A and B [7] 

2B = A[(-2 + 3n)(l + a)]. (17) 

Given A or B and n we can uniquely determine B or A, respectively. For A oc 1/t 2 , there 
is a relation between Qm and the time-dependence of scaling factor R(t). Assuming that 
VLm does not change during the matter-dominated era (a = 0) [8] 

2 

R(t) = (jy nM {tlMCtf)^. (18) 

The quantities q, R(t) and Qm are mutually related. Given one of them the remained 
quantities can be determined by means of Eqs. (18), (16) and (15). It seems that Qm is 
best determined by the observation, therefore, we take it for the calculation of q and R(t). 
Inserting Qm = 1/3 into Eq.(18) we obtain R oc t 2 which yields q = — l/2 = l/n — 1, 
n = 2. We see from Eq.(9) that B = 1 which inserting into Eq.(17) gives A = (1/2). 
The mean energy density p and the cosmological constant is given as (a = 0) 

P = T^ and A= 8^' (19) 



respectively. Their ratio 
Supposing the flat space, we have 



j = 1/2- (20) 



1 2 
Q M = - and Q\ = -■ (21) 

There is no "fine tuning" problem in our model since the ratio of the A-part energy 
density to the mass-energy density of the ordinary matter remains during the cosmic 
evolution constant. 

4 The force exerting on the moving bodies in the cosmic 
quaternionic field 

When the ordinary matter is present in the $-field then we have to apply the general 
field equations (3) and (5) The ordinary matter changes the value of the field variable 
$ in comparison to the pure <&-field. Simultaneously, <I>-field exerts a force f\ on the 
moving matter. This Lorentz-like force is given, similarly as in electrodynamics, by the 
expression 

fj = ^Fij- (22) 
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When calculating this Lorentz-like force we take, for the save of simplicity, only the 
pure $-field, i.e. that <E>-field in which no ordinary matter is present, j\ m ^ = -^Pm-ofi- 
Inserting J- m ^ in Eq. (22) we find for the Lorentz-like force acting on the moving mass 
body in the presence of the $-field the following formula [15]. 

Fi = c^VGmo^u, = c^VG^pi. (23) 

This is the fundamental formula which we apply in our further consideration. 

It is to be expected that the cosmic <I>-field manifests itself in the present-day only 
in the following astrophysical situations: (i) at the large mass concentrations, (ii) at the 
large velocities of massive objects (iii) during the large time and space scales. For the 
sake of simplicity, we confine ourselves to the non-relativistic case, i.e. we suppose that 
m = const, and v <C c. Then Eq.(6) turns out to be 

mv = c^VG&mv. (24) 

Below, we present three possible effects of the quaternionic field in the cosmic conditions: 

(i) The increase of the velocity of the moving bodies in the Q -field. Since c -1 \/G$ = 1/t 
we get a simple differential equation v = 0v where (5 = 1/t the solution of which is 
v = Ct. A free moving object in the quaternionic field is accelerated by a constant 
acceleration C. This acceleration is due to the immense smallness of (3 1/10 18 in 
the present-day extremely small. As is well-known for a given time instant the Hubble 
constant H is equal in the whole Universe. Supposing [3\q = H, the solution of equation 
v = Hv becomes v = Hr + C, where C is an integration constant. Setting C = we get 
a Hubble-like law v = Hr. 

(ii) The increase of the kinetic energy of the moving bodies in the -field. The gain of 
kinetic energy of a moving body per time unit in the quaternionic field is (fi || Vi) 

— = FiVi = c-'VGPmv 2 = 2^Gc- 1 ^E km = 2(3E kin . (25) 

Again, the increase of the kinetic energy of a moving object is extremely small. However, 
for a rapid rotating dense body it may represent a considerable value. For example, 
a pulsar rotating around its axis with the angular velocity lo having the moment of 
insertion /. Its kinetic energy is E kin rj Iuj 2 and its change in the quaternionic field is 
dEki n /dt f» filJ 1 w 10 32 which is a value only of some orders of magnitude smaller than 
the energy output of a pulsar [36]. 

(iii) The change of the kinetic parameters of the gravitationally bounded moving bodies. 
This can be best demonstrated by describing the motion of the Earth around the Sun 
taking into account Eq.(13). It holds 

F 1+ F 2 = -^p*r, (26) 
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where F\ = m®r and F2 = m^fir. Inserting F\ and F2 in Eq.(26) we have 
from which it follows 



l3r 2 4> = const. = h (27) 
The reciprocal radius u satisfies the equation 



d 2 u 2 GMq 



the general solution of which is 



C 

u(4>) = ^ + C2 cos(0) - ci sin(</>) 
Setting ci = C2 = 0, i.e. supposing that the orbit is circle we get 

r ~ Kt 2 

Hence, the distance of the Earth and the Sun varies with time like 



4 ~ t2 - ( 2 8) 



CM. i 2 f3 2 

According to Eq.(28) the distance between the Earth and the Sun is increasing direct 
proportional to the square of time. There have been several atomic time measurements 
of the period of the Moon orbiting around the Earth. A description of the work of some 
independent research groups can be found in Van Flander's article [34]. We simply point 
that after subtracting the gravitational perturbative (tidal) effects, Van Flanders gives 

|! = - = (3.2±l).l(r 10 /yr, 
where n = 2ir/P is the angular velocity. Using Eq.(ll) we get comparable value 

\h/n\ w 4.10 10 /vr- 

However, given the complexity of the data analysis, we must certainly await further 
confirmation by different, independent test before concluding whether the $-field really 
affects the motion of cosmic bodies. Nevertheless, it can be asserted that at present, 
there exists no evidence against the influence of the <E>-field on the moving bodies at the 
level of the supposed present-day intensity of the <E>-field. Due to large value of the cosmic 
time the present-day effects of the <£-field lei on the limit of the observability. However, 
they had, probably, strong influence in the early universe . For example, the strong 
3>-field can destabilize the large rotating mass concentration (e.g. quasars) forming from 
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them the present-day galaxies. We note that the enlargement of distance of the Earth 
and the Sun is also suggested by the large numbers hypothesis presented by Dirac in 
1937. In this hypotheses Dirac supposed that G oc 1/t. The astrophysical and geological 
consequences of this hypothesis are discussed in details in [35] [33]. 

It is noteworthy that the force exerting on cosmical body in cosmic quaternionic 
field is always parallel to the direction of the velocity. This means that velocity of 
moving bodies in cosmic 3>-field is in all direction increasing. For example, the moving 
bodies stemming from a exploding cosmic body is equally speedup as those falling in the 
collapsing center. 

5 The creation of particle in the <£>-field 

Another interesting effect of the cosmic <3?-field is the possibility of the creation of real 
particles from the virtual ones. Particle creation in nonstationary strong fields is well- 
known phenomenon studied intensively in seventies ( see, e.g. [32]). There are several 
proposed ways for the creation of real particles from the virtual ones in the very strong 
and nonstationary gravitation field. We propose here a new mechanism of the creation 
of real particles from the virtual ones in the presence of the <3?-field. 

We note that our further consideration on the creation of matter from the vacuum 
quantum excitation are done by a semiclassic way, although we realize that they should be 
performed in terms of an adequate theory of quantum gravity. However, as is well-known, 
when constructing a quantum theory of gravity one meets conceptual and technical 
problems. The usual concepts of field quantization cannot be simple applied to gravity 
because standard field quantization (e.g., elm field) is normally done in flat spacetime. It 
is impossible to separate the field equations and the background curved space because the 
field equations determine the curvature of spacetime. Moreover, the classical quantized 
field equations are linear and that of gravitation are non-linear and weak-field linearized 
gravitation field is not renormalizable. There is even no exact criterion on which time and 
space scales one has necessary to apply quantum laws for gravitation field. Therefore, we 
take the inequality \A\ < h, where A is the classical action, as a criterion for a possible 
application of quantum physics in gravity. 

According to quantum theory, the vacuum contains many virtual particle-anti-particle 
pairs whose lifetime At is bounded by the uncertainty relation AEAt > h The proposed 
mechanism for the particle creation in the $-field is based on the force relation (23). 
During the lifetime of the virtual particles the Lorentz-like force (23) acts on them and 
so they gain energy. To estimate this energy we use simple heuristic arguments. As is 
well-known, any virtual particle can only exist within limited lifetime and its kinetics is 
bounded to the uncertainty relation ApAx > h. Therefore, the momentum of a virtual 



9 



particle p is approximately given as p ~ hAx" 1 . If we insert this momentum into Eq.(23) 
and multiply it by Ax, then the energy of virtual particle AE, gained from the ambient 
<&-field during its lifetime, is 

FAx = AE = VG$(t)^. (29) 

When the $-field is sufficiently strong then it can supply enough energy to the virtual 
particles during their lifetime and so spontaneously create real particles from the virtual 
pairs. The energy necessary for a particle to be created is equal to m v c 2 (m v is the 
rest mass of the real particle). At least, this energy must be supplied from the ambient 
<3?-field to a virtual particle during its lifetime. Inserting <3? into Eq. (29), we have 

h 

AE ~ l 7> 

(t + to) 

Two cases may occur: (i) If m v c 2 < AE, then the energy supplied from the $-field is 
sufficient for creating real particles of mass m v and, eventually, gives them an additional 
kinetic energy, (ii) If m v > c 2 AE, then the supplied energy is not sufficient for creating 
the real particles of mass m v but only the energy excitations in vacuum. 
The additional kinetic energy of the created particles, when AE > m v c 2 , is 

E kin = AE - m v c 2 = h - m v c 2 . 

(t + 1 ) 

During the time interval (ss 0, 10~ 20 ) after the Big Bang, the masses of the created 
particles lei in the range from 10 -5 to 1CT 27 g. Their kinetic energy was E^ n = [h/ (t + 
to)] — moc 2 . reached values up to 10~ 5 erg, which corresponds to the temperature of 10 21 
K. Today, energies of the virtual pairs, gained during their lifetime, are negligible small. 



6 The energy density of the cosmic $-field as a possible 
candidate for the dark energy 

As has been shown in the previous Chapters, when taking the field energy density of 
cosmic quaternionic field as the vacuum energy density, the problems presented in the 
Introduction may be resolved. The problem of the cosmological constant, because the 
value of A is consistent with data, the tuning and age problems because the ratio of mass 
to vacuum energy density does not vary during the cosmical evolution and the age of the 
universe is large enough to evolve the globular cluster. The flatness problem is also solved 
because the sum of the dark energy to the ordinary matter yields just the critical energy 
density. Moreover, in the cosmic quaternionic field, there exists a plausible mechanism 
of the particle creation. We conclude that the energy density of the <3?-field represent 
a relativistic quantity satisfying Gliner's requirements which is smoothly distributed in 
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space. It causes the speedup of the universe, balances the total energy density to the 
critical one and gives a plausible mechanism for the particle creation. 

It is generally believed that dark energy was less important in the past and will be- 
come more important in the future. In our model the value of dark energy is proportional 
to 1/i 2 , therefore, its value becomes very large in past and will be adequate large in the 
future. We remember that the force exerted by the cosmic <I>-neld on moving bodies 
acts always in the direction of the velocity. That means that the high value of the dark 
energy in the early universe does not interfere with the structure forming, contrarily, it 
accelerates it. 

The evolution of the universe with the cosmic quaternionic field can be briefly 
sketched as follows: The cosmic evolution started purely field-dominated era with the in- 
flation, after which a massive creation of particles began together with enormous release 
of entropy. The masses of the created particles reaches values up to —5 g and the kinetic 
energy of the created nucleons values up to 10~ 5 erg, which corresponds to the temper- 
ature of 10 21 K. The large vacuum energy density of the cosmic quaternionic field at the 
early stage of the universe accelerates its structure formulation. From what has been 
said above we conclude that the energy density of the cosmic quaternionic field might 
be a possible candidate for the dark energy because (i) it has the value consistent with 
data (ii) it does not suffer from the cosmological constant, fine-tuning, age and flatness 
problems (iii) it yields a plausible mechanism particle production and (iv) it accelerated 
the structure formulation in the early universe. 

Motivated by the desire to find a possible candidate of the dark energy among the 
family of quaternionic fields we found the suprisingly simple quaternionic field whose 
energy density might be considered as the dark energy. (This points out that also the 
classical field may be interested by study of the quintessence [41]. The energy density 
of this field has the desired properties of dark energy and changes generally our view 
of the vacuum energy density modelled by the cosmological constant. It could not be 
seen as the carrier of repulsive gravity but as an amplifier of velocity of the moving 
bodies independently of the direction of their motion. The fact that the vacuum density 
accelerates the expanding galaxies is cause do to fact that they move from the center 
of the universe. When the galaxies would collapse the vacuum energy would accelerate, 
likewise, their collapsing. This may eventually lead to change the basic equation of the 
Friedmann cosmology. As well-know there are many versions of Mach's principle in the 
literature and a unique, satisfactory formulation of this principle does not seem to exist as 
yet. It seem that the influence of the moving bodies changes the value the field variable $ 
which for its part affects other cosmic bodies might be seen as a form of Mach's principle. 
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APPENDIX 



The quaternionic field equations can be described by a quaternionic equation con- 
sisting of the quaternionic differential operator, the field and source quaternions [5]. 

(i) The quaternionic differential operator is the quaternion (i,j,k are the quater- 
nionic units obeying the following relations ij = —ji, ik = —ki, jk = —kj and 
e = f = k 2 = -l) [5] [9] 

i— i -d .3 , d s d . — - 

(ii) The field quaternion the components of which are the field variables is 

* = i*i + + k^ 3 + * 4 , (A2) 

where 

*1 = * 2 = ($2 + S$ 2 ); ^3 = ($3 + S$ 3 ), ^4 = ($4 + S$ 4 ). 

(hi) The source quaternion 

47T 

J = — [t J\ + jJ 2 + kJ 3 + J 4 ] , (A3) 
c 

where 

J\ = {j' x - sJ x), J2 = {Jy-sJ y ) J 3 = (J z - sJ z )}, J4=jo + sjo- 

The general quaternionic field is described by the field equation of the type (see, e.g. 
[10]) _ 

□ d> = J (A4) 

The field quaternion $ consists of the vector part = <3?i + s$2> where $i = 

($1, $ 2 , $3), $2 = ($1, $2, $3) and the scalar part ^4 = ($4 + S&4). 

The source quaternion consists likewise of the vector part J = J\ — sJ 2 , where 

J*i = (J^., J z , j' z ) and J2 = (J x , J y , Jz) and the scalar part J4 = 47r(jo + sj' Q ). 

The energy density of a quaternionic field is [5] 

4 

E = £<M>*- (A5) 
fc=i 

By using the vector notation we can rewrite the quaternionic field equations (A4) in the 
form 

-—i + v.* = 4vrJ 4 
cot 

^--V$ 4 + Vxf = -J. (A6) 
cot c 
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According to the specification of the field variables and source components in the 
field and source quaternion, respectively, we get the following fields: 

(i) If we choose <J> = —E, = —B, J = — J j = -g, <J> 4 = <S>' 4 = f = j' = 
and we associate J, jo with the components of electromagnetic 4-current, we get just the 
standard Maxwell equations [30] [31]. 

(ii) If we associate, in addition, j' and j' with the components of the monopole 4-current 
then we get the Maxwell equations with monopoles [29]. 

(v) If we associate the components of the field quaternion with the field variables as in 
(i) and take J, jo as electric and J , j as monopole currents, $4 and <I> 4 as the scalar and 
the pseudoscalar variables, respectively, then we obtain the Ohmura field equations [18]. 

(vi) If we put $ = <P = $'4 = J' = f = 0, and take J*, j* as not specified 4-current, we 
get the scalar quaternionic field equations. 

We note that all of these field equations can be written also in form of tensor equation 
by means of the corresponding field tensors [5]. The scalar quaternionic field represents 
the $-field, whose field equations in vector form we consider in Chapter 2. 
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